Search for solar axions in XMASS, a large liquid-xenon detector 
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Abstract 

XMASS, a low-background, large liquid-xenon detector, was used to search for solar axions that 
would be produced by bremsstrahlung and Compton effects in the Sun. With an exposure of 
5.6 ton days of liquid xenon, the model-independent limit on the coupling for mass <C 1 keV is 
Idaeel < 5.4 X 10^"*^^ (90% C.L.), which is a factor of two stronger than the existing experimental 
limit. The bounds on the axion masses for the DFSZ and KSVZ axion models are 1.9 and 250 eV, 
respectively. In the mass range of 10-40 keV, this study produced the most stringent limit, which 
is better than that previously derived from astrophysical arguments regarding the Sun to date. 
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The axion is a hypothetical particle invented for solving the CP problem in strong inter- 
actions As the initial Peccei-Quinn-Weinberg-Wilczek model of axions is directly tied to 
the electroweak symmetry-breaking scale, an experimental search was relatively easy and the 
model was ruled out early. However, invisible axion models such as DFSZ j2| and KSVZ sl, 
whose symmetry-breaking scale is separated from the electroweak scale, are still viable. The 
DFSZ axions have direct couplings to leptons whereas the KSVZ axions (hadronic axions) 
do not have tree-level couplings to leptons. In these models, the mass of axions is 

a/z /^m^ 6.0 eV 

777 = = 

" l + Z fa /a/106 GeV 

where /„, /tt, and are the axion decay constant ^, the pion decay constant, and pion 
mass, respectively, and z = nid/mu ~ 0.56 is the quark mass ratio. 

At present, the search for axions as well as axion-like particles (ALPs) focuses on cou- 
plings to photons (ga-y-y), uuclcons (gaNN) and electrons (gaee)- There are three types of 
searches: (1) laboratory-based experiments in which sources and detectors are prepared, (2) 
astrophysical investigations that examine any significant deviations in the properties of stars 
from theoretical predictions due to extra emission of energy, and (3) using laboratory de- 
tectors to look for axion signals from the Sun or cosmological relics. Experiments searching 
for axions have so far produced null results, but sensitivities continue to improve. 

In experimental searches that utilize of experiments using strong magnets 



-[7| successfully improved sensitivities by increasing the magnetic field strength and the 
conversion length. The suggestion [8] to use Bragg scattering to improve sensitivity for solar 



axions in crystalline detectors was used in IJllSl. Another way to enhance sensitivity is to 
exploit resonant absorption on nuclei 13|. To date, several experimental results are obtained 



in this scheme 
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24j |. Significant improvement can be achieved if the signals can be read 



out efficiently. On the other hand, an efficient experimental search with gaee has not been 
performed. A pioneering experiment used a Ge detector (710 g) 25| and a recent search used 
a Si(Li) detector (1.3 g) to search for signals from axions generated by the bremsstrahlung 
and Compton effect via the axioelectric effect 26 1. 

The choice of target material strongly affects the reach of a solar axion experiment using 
axion coupling to electrons. Liquid Xe is both dense and has a high atomic numbers 
The XMASS detector, which uses 835 kg of liquid xenon in its sensitive volume, is suitable 
for this purpose. Its low energy threshold (0.3 keV) is also useful as the predicted energy 
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29|. The 



spectrum is very soft and has a peak at less than 1 keV for hght axions. Its low background 
(a few keV~^kg^^day^^) makes it particularly useful when searching for solar axions. 

The signals we searched for are produced by the Compton scattering of photons on 
electrons e + 7 — > e + a and the bremsstrahlung of axions from electrons e + Z— >e + a + Z 
in the Sun. The expected fluxes and spectra are derived as follows. 

The solar axion flux produced by Compton scattering was calculated in 
axion differential flux is expressed as 

dEa Jo JEa dE^ dEa ^ ^ ' ^ ' 

where Ea is the total energy of the axions, A is the average distance between the Sun and 
the Earth, Rq is the radius of the Sun, dN^/dE^ is the blackbody spectrum of photons, 
da'^ /dEa is the cross section for the Compton effect, and N(.{r) is the electron density at the 
radius r. Since and E^ is assumed to be much smaller than me, the differential cross 
section is approximately a product of 6{Ea — E^), and the total cross section 29| is expressed 
as 
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where a is the fine structure constant, mg is the electron mass, gaee is the axion's coupling 
to electrons j4| which is (1/3) (cos^ /3)me//a in the DFSZ axion model 26|, and fa = (1 — 
m'l/E^Y^'^ is the velocity of the outgoing massive axion. cot (3 is the ratio of the two Higgs 
vacuum expectation values of the model 4|. 

The energy spectrum of solar axions produced by the bremsstrahlung effect was calculated 
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301]. The differential energy spectrum is 



d^ 

dEn 



«0/-~ dN, da'' 



dEr 



dEeY,Z^N{rydr, 



(3) 



A"^ Jo JEa dEe z,A 

where is the velocity of the electrons, dNe/dE^ is the energy spectrum of the electrons, 
da'' /dE^ is the cross section for the bremsstrahlung effect, and A^z.a(^) is the atom density 
at radius r. The cross section da'' I dEg is calculated by considering the energy conservation 
of the electron and axion system 30 1. 

The temperature, electron density, and atomic density are given by the standard solar 
model BP05(OP) 31]- Figure 1 in Ref. {26] shows the energy spectra for various masses 
of axions. The bremsstrahlung component dominates below lOkeV, whereas the Compton 
contribution dominates at higher energy. 




energy (keV) 



FIG. 1. Expected energy spectra of events observed using the liquid-xenon detector. No resolution 
effects are included. Different curves are for axion masses with 0, 1, 2, 4, 8, and 16keV. The 
inset shows spectra of axion masses with 32 and 64keV. Due to a cross section enhancement for 
nonrelativistic axions, an increase at ~ rria can be seen. The step around 5 keV corresponds to 
the L-shell absorption edge of the axioelectric effect. 



The expected energy spectrum to be observed with a detector is 



dE 



dEn dEn 



(4) 



Ea = E 



where aae{Ea) is the cross section for the axioelectric effect j32l|. For the cross section, the 
expression of Eq. (3) in Ref. |26| is used for Va 
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O'aeiEa) = (JpeiEa 
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where ape{Ea) is the photoelectric cross section of the detector medium for gamma rays with 
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34| . The predicted energy 



energy Ea- The photoelectric cross section is available in Ref. 
spectra for a xenon target for various axion masses are shown in Fig. [H 

The predicted energy spectra calculated above are used to generate Monte Carlo sim- 
ulation samples. Axion signal samples can be simulated by injecting gamma rays whose 
energy is the same as the total energy of the incoming axions. This is because (1) there is a 
relationship between the cross section of the axioelectric effect and the photoelectric effect 



as in Eq. ([5]), (2) the photoelectric effect is dominant in this energy range (<100keV), and 
(3) the process after the axioelectric effect is exactly the same as that for the photoelectric 
effect. In the simulation, we considered the nonlinearity of the scintillation yield for gamma 
rays, the optical processes of the scintillation photons in the detector, the photoelectron 
distributions and discrimination threshold of photomultipliers, and the trigger conditions of 
the data acquisition system. The detailed description of the simulation and efficiencies were 



previously reported 
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36|. After taking into account the reduction efficiency described in 
the next section, the expected energy spectra for various masses of axions are obtained. 

The XMASS detector is a large liquid-xenon detector located underground (3000 m water 
equivalent) at the Kamioka Observatory, Japan. It contains an 835-kg liquid-xenon target 
with a surface of a pentakis-dodecahedron that is tiled with inward looking photomultiplier 
tubes (PMT), 630 of which have hexagonal and 12 have round photocathodes. The PMTs 
(R-10789, Hamamatsu) are specially developed for this low-background detector. The pho- 
toelectron yield at the center of the detector is evaluated at 14.7photoelectrons (p.e.)/keV 
using an internal ^^Co source. The positional dependence (maximum 15%) of the photo- 
electron yield caused by the angular acceptance of PMTs and absorption of scintillation 
light are taken into account in the Monte Carlo simulations. Data acquisition is triggered if 
four or more PMTs have more than 0.2 p. e. within 200 ns. The trigger efficiency around the 
trigger threshold was examined by LEDs placed at the detector wall. The observed behavior 
was well reproduced by the Monte Carlo simulations. Signals from each PMT are fed into 
charge ADCs and TDCs whose resolution is around 0.05 p. e. and 0.4 ns, respectively. The 
liquid-xenon detector is surrounded by a water Cherenkov veto counter, which is 10.5 m in 
height and 10 m in diameter. It is equipped with 72 20-inch PMTs whose signals are fed 
into the ADCs and TDCs. Data acquisition is triggered if eight or more 20-inch PMTs have 
hits. The detector is described in detail in Ref. 



The data set used in the solar axion search experiments covers February 21-27, 2012. 
A sequence of standard data reduction is applied to remove events caused by afterpulses 
and electronic ringing. The standard reduction consists of a series of cuts: (1) the event is 
triggered only by the liquid-xenon detector; (2) the time difference to the previous event is 
more than 10 ms; (3) the root mean square of the hit timing is less than 100 ns and is used to 
reject events caused by afterpulses of PMTs due to bright events; and (4) the number of PMT 
hits in the first 20 ns divided by the total number of hits is less than 0.6 for events in which the 
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scaled energy (keV) 



FIG. 2. Observed energy spectra. The horizontal axis shows the "scaled energy" calculated by 
dividing the number of photoelectrons by the photoelectron yield at the center of the detector, 
14.7p.e./keV. Error bars are statistical only. In the same figures, the efficiencies (open circles with 
horizontal bars; 1 for 100% efficiency) for the Cherenkov cut (see text) are also shown. The inset 
shows those at high energies. 



number of photoelectrons is less than 200. The fourth cut was applied to remove Cherenkov 
events originated from in photocathodes (Cherenkov cut). The energy threshold of this 
analysis is low (0.3 keV) because of our exceptional photoelectron yield, which is the largest 
among current low-background detectors. A more detailed description of the reduction can 



be found in Ref. 



35|. 



Figure [2] shows the observed energy spectra. The total livetime is 6.7 days after consider- 
ing the dead time caused by the cut (2). The effect of trigger cut (1) is visible below 0.4 keV 



as shown in Fig. 3 in Ref. 35| and is considered in our Monte Carlo simulations. The same 
samples show that the cut (3) has negligible effect on the signals. The signal efficiency due 
to the Cherenkov cut, which is drawn in the same figure, was conservatively evaluated using 
low-energy gamma-ray sources such as ^^Fe and ^^^Am sources at various positions. Because 
the efficiency weakly depends on the radial position of the events and gradually decreases 
outward, the efficiency adopted in the analysis was mostly evaluated at a radius of 40 cm 




FIG. 3. Comparison between the observed data (points with error bars) and expected spectrum 
(sohd histogram) for axion masses of 0, 5, 10, and 50keV. The sohd histograms are scaled to the 
maximum couphng ahowed at 90% C.L. 

where 93% of the mass was contained inside. The Monte Carlo samples were compared with 
the observed energy spectra after weighting this efficiency. 

The observed spectra do not have any prominent features to identify axion signals with 
respect to the background. Instead, strong constraints on Qaee can be obtained from the 
observed event rate in the relevant energy range. In order to set a conservative upper limit 
on the axion-electron coupling constant (yf^ee, the coupling is adjusted until the expected 
event rate in XMASS does not exceed the one observed in any energy bin above 0.3 keV. 
Figure |3] shows the expected energy spectra with the coupling constants obtained by the 
procedure above. Figure H] shows the summary of the bounds of Qaee- For small axion 
masses, a Qaee value of 5.4 x 10^^^ is obtained. This is the best direct experimental limit to 
date and is close to that derived from astrophysical considerations based on measured solar 



neutrino fluxes: Qaee = 2.8 x 10~^^ 38|]. For axion masses > lOkeV the energetics in the 
Sun are no longer sufficient to effectively produce such axions. A systematic uncertainty 
inherent to our method of comparing bin contents arises from the specific choice of binning. 
This and systematic uncertainties for energy scale including energy threshold, Cherenkov 
cut efficiency, and energy resolution are evaluated to be 2%, 1%, 2%, and 1%, respectively. 
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FIG. 4. Limits on gaee- The thick sohd hne shows the hmit obtained in this study. The other sohd 

25| . Si(Li), ^^^Tm, reactors, o-Ps, beam- 



lines are hmits obtained by laboratory experiments: Ge 
dump experiments, dark matter search experiments. The last ones assume dark matter consists of 
axions. (See [2Cj_33] and references therein.) The dash-dotted lines show astrophysical limits from 



red giant stars [4] and the solar neutrino flux 



The dashed lines are theoretical predictions for 



the DFSZ (cos^ /? = 1) and KSVZ {E/N = 8/3) models. This study gives a stronger constraint 
by a factor of two over previous direct experimental limits for axion mass ^IkeV, and the best 
constraint absolute between 10 and 40 keV. 



The total systematic error, 3%, is obtained by summing these contributions in quadrature, 
and the limit in Fig. H] (90% CL.) takes this error into account. 

The calculated limit depends on the interaction processes considered in our detector as 
well as the processes considered for solar axion production in the Sun. Processes such as the 
inverse Primakoff effect and nuclear absorption on the detection side, and the Primakoff effect 
and nuclear deexcitation on the production side can be neglected because the constraints on 
gg^j and gaNN are tight. A possible additional contribution caused by gaee on the detection 
side is the inverse Compton effect. This can be neglected because of its small cross section 



39| . On the production side, there are other known contributions such as electron-electron 



bremsstrahlung 



40j and the axio-recombination effect 



4l|. However, the expected fluxes 



for these processes are only known in the hmit of massless axions. For this reason and in 
order to directly compare our results with the most relevant previously published ones we 
restrict the production processes we consider to the electron-nuclei bremsstrahlung and the 
Compton effect. As omitting production mechanisms lowers the flux estimate, all the limits 
thus derived will have to be considered conservative. 

The nature of the events surviving the analysis cuts is also of interest. According to 



our study on these events, most of them originate on the inner surface of the detector [42 1. 
These events are attributed to radioactive contamination in the aluminum seal of the PMT 
entrance windows, ^^C decays in the GORE- TEX® sheets between the PMTs and the 
copper support structure, and light leaking from gaps in between the triangular elements of 
this support structure. 

In summary, solar axions produced through axion-electron coupling were searched for 
in XMASS, a large liquid-xenon detector. The energy threshold is low (0.3 keV) because 
of our exceptional photoelectron yield, which is the largest among current low-background 
detectors. As our observed spectrum does not show any indications of axion signals, we 
derive constraints on the Qaee coupling. Our limit on gaee for axions with mass much smaller 
than IkeV is 5.4 x 10^^^. The bounds on the axion masses for the DFSZ and KSVZ axion 
models are 1.9 and 250 eV, respectively. For axion masses between 10 and 40keV, our new 
limits are the most stringent that are currently available. 

We gratefully acknowledge the cooperation of Kamioka Mining and Smelting Company. 
This work was supported by the Japanese Ministry of Education, Culture, Sports, Science 
and Technology, Grant-in- Aid for Scientific Research, and partially by the National Research 
Foundation of Korea Grant funded by the Korean Government (NRF-2011-220-C00006). We 
thank J. Redondo for useful discussion. 



* Current address: Research Center for Neutrino Science, Tohoku University, Sendai, 9800-8578, 
Japan 

[1] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett 38, 1440 (1977); R. D. Peccei and H. R. Quinn, 
Phys. Rev. D 16, 1791 (1977); S. Weinberg, Phys. Rev. Lett. 40, 223 (1978); F. Wilczek, Phys. 
Rev. Lett. 40, 279 (1978). 

10 



[2] M. Dine, W. Fischler, M. Srednicki, Phys. Lett B 104, 199 (1981); A. R. Zhitnitsky, Sov. J. 

Nucl. Phys. 31, 260 (1980) [Yad. Fiz. 31, 497 (1980)]. 
[3] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, 

Nucl. Phys. B 166, 493 (1980). 
[4] G. G. Raffelt, Lect. Notes Phys. 741, 51 (2008). 
[5] D. M. Lazarus et al, Phys. Rev. Lett. 69, 2333 (1992). 

[6] S. Moriyama et al, Phys. Lett B 434, 147 (1998); Y. Inoue et al, Phys. Lett. B 536, 18 
(2002); Y. Lioue et al, Phys. Lett. B 668, 93 (2008). 

[7] S. Andriamonje et al, JCAP 0704, 010 (2007); E. Arik et al, JCAP 0902, 008 (2009); 
S. Aune et al, Phys. Rev. Lett. 107, 261302 (2011). 

[8] E. A. Paschos and K, Zioutas, Phys. Lett. B 323, 367 (1994). 

[9] F. T. Avignone et al, Phys. Rev. Lett. 81, 5068 (1998). 
[10] R. Bernabei et al, Phys. Lett B 515, 6 (2001). 
[11] A. Morales et al, Astropart. Phys. 16, 325 (2002). 
[12] Z. Ahmed et al, Phys. Rev. Lett 103, 141802 (2009). 
[13] S. Moriyama, Phys. Rev. Lett. 75, 3222 (1995). 
[14] M. Krcmar et al, Phys. Lett. B 442, 38 (1998). 

[15] M. Krcmar, Z. Krecak, A. Ljubicic, M. Stipcevic, and D. A. Bradley, Phys, Rev. D 64, 115016 
(2001). 

[16] K. Jakovcic, Z. Krecak, M. Krcmar, A. Ljubicic, Rad. Phys. and Chem., 71, 793 (2004). 

[17] A. V. Derbin et al, JETP Lett. 85, 12 (2007). 

[18] T. Namba, Phys. Lett. B 645, 398 (2007). 

[19] A. V. Derbin et al, Buh. Rus. Acad. Sci. Phys. 71, 832 (2007). 

[20] P. Belh et al, Nucl. Phys. A, 806, 388 (2008). 

[21] A. V. Derbin et al, Eur. Phys. J. C 62, 755 (2009). 

[22] A. V. Derbin et al, Phys. Atom. Nuclei 74, 596 (2011). 

[23] G. Bellini et al, Phys. Rev. D 85, 092003 (2012). 

[24] P. BelH et al, Phys. Lett. B 711, 41 (2012). 

[25] F. T. Avignone, R. L. Brodzinski, S. Dimopoulos, G. D. Starkman, A. K. Drukier, 

D. N. Spergel, G. Gelmini, B. W. Lynn, Phys. Rev. D. 35, 2752 (1987). 
[26] A. V. Derbin et al, JETP Lett. 95, 379 (2012), arXiv:1206.4142| y2. 

11 



[27] F. T. Avignone III, R. J. Creswick, S. Nussinov, Phys. Lett. B 681, 122 (2009). 

[28] A. V. Derbin, A. S. Kayunov, V. V. Muratova, D. A. Semenov, and E. V. Unzhakov, Phys. 

Rev. D 83, 023505 (2011). 
[29] M. Pospelov, A. Ritz, and M. Voloshin, Phys. Rev. D 78, 115012 (2008). 
[30] A. R. Zhitnitskh and Y. I. Skovpen', Sov. J. Nucl. Phys. 29, 513 (1979). 
[31] J. N. Bahcah, A. M. Serelh, and S. Babu, Astrophys. J. 621, L85 (2005) 
[32] A. Derevianko, V. A. Dzuba, V. V. Flambaum, M. Pospelov, Phys. Rev. D 82, 065006 (2010). 
[33] M. J. Berger, J. H. Hubbel, S. M. Seltzer, et al., XCOM: Photon Cross Sections Database 

http : // www. nist . gov /pml/dat a/ xcom / index, cfm 
[34] WM. J. Veigele, Atomic Data Table 5, 51 (1973). 
[35] K. Abe et al. (XMASS collaboration), Phys. Lett B 719, 78 (2013). 
[36] K. Abe et al. (XMASS collaboration) . larXiv: 1301.28151 
[37] C. E. Aalseth et al, Phys. Rev. Lett. 106, 131301 (2011). 
[38] P. Gondolo and G. G. Raffelt, Phys. Rev. D 79, 107301 (2009). 
[39] R. Bernabei et al, Int. J. Mod. Phys. A 21, 1445 (2006). 
[40] G. G. Raffelt, Phys. Rev. D 33, 897 (1986). 

[41] S. Dimopoulos, J. Frieman, B. W. Lynn, and G. D. Starkman, Phys. Lett. B 179, 223 (1986). 
[42] Y. Suzuki for the XMASS collaboration, Proceedings of the International Workshop on the 
Identification of Dark Matter, Chicago, USA, 23-27 July, 2012. 



12 



